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Abstract: We compute the complete imaginary part of the NRQCD Lagrangian at order 
in the heavy-quark mass expansion, which includes center of mass operators, and 
at order in the matching coefficients. We also compute the imaginary part of the 
NRQCD Lagrangian at order and at order that contributes to the S-wave and 

P-wave inclusive decay widths of heavy quarkonium into light hadrons at order in 
the heavy-quark velocity expansion. If we count as(-^) ~ the calculation provides 
the complete next-to-leading order corrections to the P-wave hadronic widths, and in the 
original NRQCD power counting, the complete next-to-leading order corrections to the 
vector S-wave widths, and part of the next-to-next-to leading order corrections to the 
pseudoscalar S-wave widths. In the S-wave case, we confirm previous findings and add 
new terms in a more conservative power counting. In the P-wave case, our results are 
in disagreement with previous ones. Constraints induced by Poincare invariance on the 
NRQCD four-fermion sector are studied for the first time and provide an additional check 
of the calculation. Perspectives for phenomenological applications are discussed. 
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1. Introduction 

Non-relativistic effective field theories (NR EFT) of QCD ||, || hke non-relativistic QCD 
(NRQCD) 1^, ^ offer a systematic framework to access heavy-quarkonium properties and 
in particular inclusive decay widths. Decay width formulas may be organized in a double 
expansion in the strong coupling constant as , calculated at a large scale of the order of the 
heavy-quark mass M, and in the heavy-quark velocity v. Both expansion parameters are 
relatively small. In the bottomonium system, typical reference values are Qs(Mf,) ~ 0.2, 
~ 0.1 and in the charmonium one, as(Mc) ~ 0.35, ~ 0.3. 
The increasing accuracy of the experimental measurements |jl], ^, |7| calls for a corre- 
sponding accuracy in the theoretical predictions. The inclusive decay widths of J/V'; ^(25*) 
and T(IS') into light hadrons are presently known within a few percent uncertainty, while 
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the uncertainties in the inclusive decay widths of T{2S) and T(3S') are less than 10% j5|. 
Theoretical accuracies of about 5% both in the charmonium and in the bottomonium case 
require at least the calculation of 0{v^, Os ) corrections. The S-wave decay of r]c into 
light hadrons is presently known within a 15% uncertainty, while the P-wave decays of Xcj, 
with J = 0, 1, 2, are known within a 10% uncertainty [|[. In the P-wave case, the improve- 
ment of the experimental accuracy has been noticeable over the last few years and the data 
are now clearly sensitive to next-to-leading (NLO) corrections |^, |^. Hence, for the decay 
of the charmonium P-wave states, theoretical accuracies matching the experimental ones 
require the calculation of 0(f ^, as) corrections. 

In this work, we consider relativistic corrections of order v'^ and to inclusive decays 
of P- and S-wave quarkonium into light hadrons respectively. The leading-order S-wave 
decay width is proportional to the square of the wave-function in the origin and is therefore 
of order v^. The leading order P-wave decay width is proportional to the square of the 
derivative of the wave- function in the origin, and is therefore of order v^. Then, corrections 
of order to S-wave decays and of order v"^ to P-wave decays provide in both cases decay 
widths at order v'^ in the relativistic expansion. We consider only processes where the 
quark and antiquark annihilate into two gluons. Hence, more precisely, the paper provides 
the Ogt''' terms of the S-wave and P-wave inclusive decay widths. 

In the S-wave case, corrections of order v"^ and were first considered in [Q] and 
Q respectively. We agree with their results if we use their power counting, but find 
additional contributions in the more conservative counting that we adopt. In the P-wave 
case, corrections of order v"^ were first calculated in jl^. Our results disagree with those 
results. In particular, we find different matching coefficients for the dimension 10 operators. 
Moreover, also adopting the power counting of our decay widths appear to contain 
two matrix elements more. 

The paper is organized in the following way. In section |2|, we set up the formalism, 
discuss the power counting, introduce our basis of operators and give the general form of the 
decay widths at order v'^ . In section |3|, we calculate the short-distance imaginary parts of 
the NRQCD four-fermion operators by matching annihilation diagrams of order ctg. Octet 
operators are calculated by matching diagrams with an external gluon. In section ^, we 
show how Poincare invariance is realized in the EFT in the form of exact relations among 
matching coefficients. Such relations provide an additional and independent check of some 
of the results. In section |5|, we conclude by summarizing the present knowledge about 
inclusive decays and discussing phenomenological applications and future developments of 
this work. In appendix ^ and we list all the operators and the matching coefficients 
that have been employed through the paper. 

2. Hadronic decay widths in NRQCD 
2.1 NRQCD 

The main mechanism for quarkonium to decay into light hadrons is quark-antiquark an- 
nihilation. It takes place at a scale which is twice the heavy-quark mass M. Since this 
scale is perturbative, quark-antiquark annihilation may be described within an expansion 
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in the strong coupling constant Og- Experimentally, this is shown by the narrow widths of 
quarkonia below the open flavor threshold. The bound state dynamics, instead, is charac- 
terized by physical scales smaller than M, such that a perturbative expansion in as may 
not be allowed. It is however possible to take advantage of the non-relativistic nature of the 
bound state and expand in the relative heavy-quark velocity v. In an EFT language, once 
the scale M has been integrated out, the information on the decays is carried by contact 
terms (four-fermion operators) whose matching coefficients develop an imaginary part Q. 
In NRQCD, the decay widths factorize in a high-energy contribution, encoded in the imagi- 
nary part of the four-fermion matching coefficients, and a low-energy contribution, encoded 
in the matrix elements of the four-fermion operators evaluated on the heavy-quarkonium 
states. The NRQCD factorization formula for the inclusive decay width of a quarkonium 
state H into light hadrons (l.h.) is 

r(i7^1.h.) = 2^^(i/|0S|F). (2.1) 

n 

\H) is a mass dimension —3/2 normalized eigenstate of the NRQCD Hamiltonian with 
the quantum numbers of the quarkonium state H. The coefficients c^"^ can be calculated 
in perturbation theory by matching Green functions or physical amplitudes in QCD and 

(n) 

NRQCD. 0^:7 stands for a generic four-fermion operator of dimension dn, whose general 
form is ■ ■ )xxH' ' ' V' being the Pauli spinor that annihilates a quark and x the 
one that creates an antiquark. The operators (• • • ) may transform as singlets or octets 
under color SU(3) gauge transformations. In the first case, we denote the operator with 
the subscript 1, in the second with the subscript 8. A list of relevant four-fermion operators 
is provided in appendix 

It is the purpose of this work to calculate the order contributions to the c^") co- 
efficients that multiply matrix elements up to order v"^ . These involve operators up to 
dimension 10. 

2.2 Power Counting 

In the factorization formula (|2.l| ), the matching coefficients c*-"-* are series in as while the 
matrix elements {H\0^^^\H) are series in v and are, in general, non-perturbative objects. 
In NRQCD, it is not possible to attribute a definite power counting to the matrix elements 
because of the several contributing energy scales. These are the relative momentum Mv, 
the binding energy Mv"^, and the typical hadronic scale Aqc^; additional scales may enter 



at higher orders in the calculation |11]. Whatever power counting one assumes, as long as 
u <C 1, matrix elements of operators of higher dimensionality are suppressed by powers of 

V. 

The NRQCD Lagrangian is constructed as an expansion in 1/M and hence it is inde- 
pendent of the power counting. We shall adopt a power counting, however, when assessing 
the size of the different matrix elements contributing to the decay widths. We will assume 
Mv of the same order as Aqcd and adopt the following rules. Matrix elements of the 
type {H'\0\H) , where 0\H) and \H') have the same quantum numbers and color trans- 
formation properties in the dominant Fock state, scale (at leading order) like {Mv)'^~^ , d 
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being the dimension of the operator O. If 0\H) and \H') do not have the dominant Fock 
state with the same quantum numbers, then the matrix element singles out a component 
of the quarkonium Fock state that is suppressed. The amount of suppression depends on 
the power counting and on the quantum numbers. As detailed in [12|, the power counting 
we adopt implies that the octet components with quantum numbers S and L ziz 1, S and 
L, S ± 1 and L of a quarkonium state are suppressed by v with respect to the singlet 
component with quantum numbers S and L, while the components with S*, L it 2 or S" it 1, 
L lb 1 are suppressed by v'^ . 

A different counting, which seems suitable for the situation Mv'^ ~ Aqcd has been 
defined in Q and used in |^, |1C]. Since our power counting is more conservative than 
the one in Q, one may recover the expressions in the latter simply by eliminating matrix 
elements that, in that counting, would be smaller than v"^: no new matrix element or 
matching coefficient is needed. 

For a critical review and a discussion on the different power countings we refer to Q 
and references therein. 



2.3 Four-fermion operators 

The four-fermion sector of the NRQCD Lagrangian contains all four-fermion operators 
invariant under gauge transformations, rotations, translations, charge conjugation, parity 
and time inversion. They may be classified according to their dimensionality and color 
content. The analysis of the four-fermion operators involved in the hadronic decay widths 
at order v"^ closely parallels the one performed for electromagnetic decays in In the 
following, we focus on the main differences, that are mostly related to the contributions of 
color octet operators to the hadronic decay widths. The presence of color octet operators, 
acting on sub-leading components of the heavy-quarkonium Fock state, is one important 
and well known characteristics of NRQCD Q. 

We organize the four-fermion sector of the NRQCD Lagrangian according to the mass 



dimension and the color structures of the operators. In section 2.3.2, we show how the 



number of independent color singlet and octet operators may be reduced by using suit- 



able field redefinitions. In section |2.3.3| , we introduce operators proportional to the total 
momentum of the heavy quark-antiquark pair: at variance with the electromagnetic case, 
such operators contribute to the decay widths at order v'^ . In appendix |A[ we give some 
details on the construction of octet operators of higher dimension and the explicit list of all 
four-fermion operators that need to be considered at the order of accuracy we are working. 



Finally, in section 2.3.4 we use the NRQCD power counting of section 2.2 to assess the 



importance of the different matrix elements and in section 2.4 we write the general form 
of the hadronic decay widths accurate up to order v'^ . 

2.3.1 Operators from dimension 6 to dimension 10 

For dimensional reasons, four-fermion operators of mass dimension 6 can only contain 
four-fermion fields, without any covariant derivative or gluon field. The only allowed color 
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structures are Ic (8) Ic and (81 1". The color octet operator 

V'^i'xX^i"^ (2.2) 

has non vanishing matrix element between the states {{QQ)sg\ ■ ■ ■\{QQ)^9) ^ which are 
subleading components of the heavy-quarkonium Fock state. Color octet matrix elements 
are particularly relevant for P-wave decays, where they contribute at leading order in the 
power counting. 

Parity conservation forbids four-fermion operators of mass dimension 7. Four-fermion 
operators of dimension 8 can be built with two covariant derivatives or with a chromomag- 
netic field. For operators built with two derivatives, the possible color structures are Ic ® Ic 
and t"" ^t"". The construction of color singlet operators is straightforward, while some care 
has to be taken in the color octet case, because of the non-Abelian nature of the gauge 
group, see appendix The covariant derivatives involved can be proportional either to 
the relative momentum of the quark and antiquark pair, for example in an operator like 

iP^'Dx-X^'Dij, (2.3) 

or to the total momentum of the pair, like in 

v(V^x) • v(xV). 

Also, operators containing both kind of derivatives can be built, like 

V'^ (-^'^) x^X-v(xV) +H.C.. 

Operators containing the chromomagnetic field can appear with the different color 
structures (g) Ic, Ic (8) t", ® and (P^H'' t^: 

'^l^^gB ■ <?xxV + H.c. , 

rV5^"-<?t\x^i^V' + H.c., 
(^^"11)^93" ■ Bt\x^filj + H.c. . 
Operators of dimension 9 can involve a covariant derivative and a chromoelectric field, 

i^hxH^ ■9E + gE-'D )^ + Yl.c., (2.5) 



and again we have to consider all the possible color structures, as in Eq. (|2.4j ). Finally, 
dimension 10 operators may involve four covariant derivatives or two covariant derivatives 
and a chromomagnetic field or two gluon fields. To clarify our terminology, we call "singlet 
operators" the ones in which both the ingoing and the outgoing QQ pairs are singlets, as 
in (p. 3D, although any covariant derivative also contains an octet part, "octet operators" 
the ones in which both the ingoing and the outgoing QQ pairs are octets, as in ( p.2D or in 
the third and fourth lines of Eq. ( |2.4D and "singlet-octet transition operators" the ones in 
which one of the two pairs is an octet and the other is a singlet, as the first two operators 
of Eq. ( ^.4| ) or the one in Eq. (p.5|). For details on the four-fermion operator definition 
and construction see appendix 
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2.3.2 Field redefinitions 



The four-fermion basis built with all possible operators allowed by rotational and transla- 
tional invariance, gauge invariance and invariance under the discrete symmetries of QCD 
is redundant since the number of four-fermion operators may be reduced by suitable field 



redefinitions. The analysis performed in Ijl 
Through the field redefinitions 



can be extended to hadronic singlet operators. 



a 



--D 
2 



D 



.XX' 



X 



(2.6) 



it is possible, for a suitable choice of the free parameter a, to trade the operator 7i_8(^5o, 
^Pi), defined in Eq. ( |A.25D , for the linear combination of Q'i(^5o) — Qi(^5'o), defined in 
Eq. ( [A. 26 ), while, through 



X 



X 



M5 

a 
M5 



D- 



D 



,XX^ 



(2.7) 



where 



^ijlk ~ 


gijglk 

3 ' 




t(i) - 

^ijlk — 


^ijn^kln 




2 ' 




t(2) _ 
^ijlk ~ 
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3 



(2.8) 
(2.9) 

(2.10) 

^ o 

the operators T^!l{^Si,^ P), with i = 0, 1, 2, can be eliminated by a suitable choice of a and 
by redefining the matching coefficients of QU^-^i), Q'd^Si), Q[{^Si,^ Di) and Q:{{^Si,^ Di) 
(see Eqs. ( [A. 25 ) and ( |A.2(^ ) for the definition of these operators). As it was noted in 
p^ , these field redefinitions do not change the sums of the coefficients h\{}So) + h"(^So), 
h[{^Si) + Ki^Si) and h[{^Si,^Di) + h'H^Si,^ Di). 

It is also possible to exploit field redefinitions to reduce the number of octet operators. 
Consider the field redefinitions 



x^x 



M5 

a 



D XXX 



-5" 



f'x 



(2.11) 



where the definition of D^t"'x is given in Eq. ( A.13| ). Eq. ( 2.11 ) induces the following 
transformation 

a 1 



M5 N, 



■'^l-8(^-Pl/ 'S'o) 
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rVs-si-^So,^ Pi) + 7777^-^8 (^5*0), 



where A'^c = 3 is the number of colors and the operators VssCSq,^ Pi) and J^sCSq) are 
defined in Eq. ( WM )- The same field redefinitions induce the following transformation on 
the kinetic term 

(2.13) 



X 



f)2 



+ ^WSo) 



Q'^eso)) 



where in the right-hand side we have neglected operators proportional to the center of 
mass momentum of the quark-antiquark pair. Equations ( 2.12 ) and ( ^.13 ) show that the 
operators 7i_8(^-Pi/ 5*0) and Qg(^So) — Qg(^5o) are not independent and that it is possible, 
for a suitable choice of the parameter a, to trade the one for a redefinition of the matching 
coefficient of the other and of Vs-sCSq,^ Pi) and TsCSq). 

With a closely related argument, introducing the field redefinitions 



X 



X 



l<-> 



'Dn,t''xx 



(2.14) 



with T-j',], given in Eqs. 



" ijlk 



(|2^)-(^l0D and DH" defined according to Eq. (|Al^ ), it 
is possible to set the parameter a in such a way that the minimal basis of operators 
either contains the three operators Ti.^i^Pj,^ Si) or, with a different choice of a, the three 
operators 1/2{Q!^{^ Si)- Ql{^Si)), 1/2{Q!^{^ Si,"" Di) - Q'^^Su'Di)) and %y'('Si,^P) 
defined in Eqs. ( |A.25| ) and ( |A.26| ). The first set of operators is more useful in dealing with 
P-wave decay widths and we will use it in the rest of the paper. 



Note that the operator Tg}^'{^Si,'^ P) as well as the operators 7^*:g ('^S'l,"^ P) previously 
introduced and T^^g' {^Si,^ P), which is required by the matching, annihilate (create) a 
singlet QQ pair with orbital angular momentum L = 1 but with no definite value of J. 
So, in our notation, we denote the annihilated pair just with its spin and orbital angular 
momentum quantum numbers, omitting the subscript J. 

2.3.3 Operators proportional to the total momentum of the quark-antiquark 
pair 

The description of the hadronic decay widths up to order v"^ requires the inclusion of 
operators proportional to the total momentum of the quark-antiquark pair into the meson. 
By parity conservation these operators must contain at least two derivatives, so they have 
at least mass dimension 8. The two derivatives can act on the QQ pair, like in 



Viae 



(2.15) 



Since the QQ pair is a color singlet, V is an ordinary derivative. If the QQ pair is a color 
octet, we can build an operator analogous to (|t|) 



(2.16) 
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where Dab is a covariant derivative in the adjoint representation. 

Also operators containing a total derivative V and a derivative D , proportional to the 
relative momentum of the pair, can be built. In this case, since under charge conjugation 
^(V'^x) ~^ ^(V'^x) ™d tp"^ Dx ^ Dx, the operators must contain a Pauli matrix in 
order to be charge conjugation invariant. An example is the operator 

Oicm = V'^ [-'2^) X • V (xV) + H.c. . (2.17) 

As explained in section the matching coefficients of the operators of mass dimension 
8 proportional to the total momentum of the QQ pair are completely determined by the 
coefficients of the dimension 6 operators. These relations are a manifestation of the Poincare 
invariance of the effective field theory. 

2.3.4 Power counting of the four-fermion operators 



From the rules given in section 2.2, it follows that 



W^'Lj)\j^^ 0,e'^^Lj)\H{^'^^Lj)) Mv'-\ (2.18) 

where \H(^^^'^Lj)) stands for a quarkonium state whose dominant Fock-space component 
is a QQ pair with quantum numbers 5, L and J, C'i(^'^^^Lj) is a singlet four-fermion 
operator that acts on the QQ pair with spin S, orbital angular momentum L and total 
angular momentum J and d is the dimension of the operator. 

The scaling of color octet matrix elements is affected by the suppression of the Fock 
state component they act on. For example, the power counting given in section |2.2| implies 



{H{^P,)\^0,{'S,)\HeP,)) ^ Mv\ 

^ (2.19) 
{H{^S^)\j^OseSi)\H{^S,)) ^ Mv\ 

In the power counting that we adopt, the gluon field and the derivative that belong 
to a covariant derivative have the same scaling. If the gluon field selects a component of 
the quarkonium Fock state, which is suppressed, like in {H(^ Pq)\Oi(^ Pq)\H(^ Pq)) , then 
its contribution to the matrix element is subleading. If, however, the gluon field selects a 
component whose projection on the operator is not suppressed or the gluon is reabsorbed 
by other gluons in the operator, then it may happen that the gluon part in the covariant 
derivative gives to the matrix element a contribution that is larger than the one provided by 
the derivative part. For example, due to the gluons in the covariant derivatives, dimension 
10 octet operators like Vs(}Pi), Qg(^S'o) and 28(^-^2), as well as the singlet operator 
Si(^-D2), contribute to the decay width of the quarkonium state H{^So) at order v'^ . 
Similar operators contribute at order f also to the decay width of the quarkonium states 
H{^Si) and H{^Pj). 

Concerning the scaling of the singlet-octet matrix elements, in the power counting of 
section |2.2| both the chromoelectric and chromomagnetic fields scale as their mass dimen- 
sion, (Mv)'^, so the scaling of a matrix element is Mv'^~^v^, where takes into account 
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the suppression of the Fock state the operator acts on. For example, consider the matrix 
elements of the dimension 8 operators defined in Eq. ( [A. 24 ): 



{HeSo)\^SM'So,' Si)\HeSo)), (2.20) 

and 

{Hes,)\—^sM's^,' Somes,)). (2.21) 

The operator Si.g,{^ Sq,^ Si) destroys a singlet QQ pair with quantum numbers ^Sq and 
creates an octet QQ pair with quantum numbers ^Si and a gluon (and viceversa), the 
operator <Si_8(^5'i,^ 5*0) destroys a singlet QQ pair with quantum numbers ^Si and creates 
an octet QQ pair with quantum numbers ^Sq and a gluon (and viceversa). Hence, both 
matrix elements scale like Mv^. 

Equations ( |A.25 ) define octet operators of dimension 9, and since the octet Fock-space 



component is suppressed hy v, we have 



and 



{Hesi)\^ri}^'es^p)\Hesi)) ~ mv\ (2.22) 

{HePj)\^ TM'Pj,'' Si)\HePj)) ~ Mv'. (2.23) 

For the reasons discussed above, in our power counting, matrix elements of octet operators 
of dimension 9, like X'8-8(^'S'o/ -Pi)) are not necessarily negligible at order -y^ because of the 
gluons in the covariant derivatives, which may couple to other gluons in the operator and 
in the quarkonium Fock state. For instance, we have 

{HeSo)\^Vs.8eSo,' Pi)\HeSo)) ~ Mv'. (2.24) 

Matrix elements of the operator Ts(^Sq) are smaller than u because of the suppression 
induced by the Gauss law. Note that also the matrix element of the following dimension 
10 operator is negligible at order v'^: 

(F(3po)|^ i^^B ■ "Dxx^^ ■ ami'Po)) ~ Mv^ (2.25) 

Finally, we discuss the scaling of matrix elements of operators proportional to the total 
momentum of the QQ pair. We work in the frame in which the heavy quarkonium is at rest. 
In this frame, operators proportional to the total momentum of the pair have non vanishing 
matrix elements only between subleading components of the heavy-quarkonium Fock state, 
containing at least one gluon. We assume this gluon to be soft {q^,q) ~ {Mv,Mv); the 
emission of a soft gluon leaves the QQ pair with a total momentum of order Mv. Hence 
the scaling of the operators V and Dab acting on the QQ pair is ~ Mv. Consider, for 
example, the matrix element of the operator Ogcm between ^Si states 

{Hesi)\oscn.\Hesi)) = esii^h^ • (^^*'^) + +■■■■ 

(2.26) 
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The leading order contribution to the l.h.s. of Eq. (|2.26 ) comes from the matrix element 



between the components \^Si) and |("^5o)85) of \H{'^Si)), the gluon in the incoming state 
being annihilated by the gluon field in D . The matrix element in Eq. ( 2.26[ ) gets a v 



suppression from each derivative, and a further v suppression from the 1(^5*0)85) state. 
Therefore it scales like and is suppressed by with respect to the leading contribution 
to the decay width. The operator Vsacm has nonvanishing matrix element if both the 
incoming and outgoing states contain a gluon. For example, it contributes to the decay 
width ofHi^Po): 

{HePo)\Vsacm\HePo)) = ieSMDl, (v^t^Vx) DI, . (x^tW) + H.c.|(35i)85) + • • • • 

(2.27) 

The matrix element in Eq. ( |2.27| ) gets two powers of v from the derivatives and two from 
the states, so it scales like f^, and contributes to the P-wave decay width at the order we 
are interested in. 

We note that for electromagnetic decays, operators proportional to the total momen- 
tum of the QQ pair do not contribute to decay widths calculated in the quarkonium center 
of mass rest frame. The reason is the following. Electromagnetic operators are obtained by 
inserting the vacuum projector |0)(0| in hadronic operators. As a consequence, any matrix 
element involving derivatives acting on both the quark-antiquark fields may be reduced 
by integration by parts either to a matrix element that does not involve an operator with 
derivatives acting on the quark-antiquark fields or to a global derivative of a matrix element 
of the type (0|(. . . )\H). The first one is a standard matrix element that does not involve 
the center of mass momentum, the last one vanishes in the quarkonium center of mass rest 
frame. 

2.4 Hadronic decay widths 

Having assumed a power counting and having chosen a basis of operators, we are in the 
position to provide explicit factorization formulas for S-wave and P-wave inclusive decays. 
The S-wave decay widths at order f ^ are: 



^-^^^^{HeSo)\Vsacm\HeSo)) + ^^^^ {H So)\0^ ^^\H So)) 
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+^^^^^^^0^(i^(^5o)|I?8-8(^5o,iPi)|i/(i5o)), (2.28) 



r(^5i - l.h.) = ^-^^^^^^{HeS^)\0,eS,)\HeS,)) 

J=0 



J=0 



2 



J=0 



M5 
fc=0,2 

2Im58eP2,^F2)^^(35^)|^^(3p^3^^)|^(3^^)^_ (2.29) 



M6 

In Eqs. ( 2.2§| ) and ( |2.29| ), the first matrix element scales like v^, the following four in the 
second and third line like v^, the following two like and the others like f ''. S-wave decay 
widths at order v'^ were computed in p. For r(^S'o —>■ l-h.), the decay width in does not 
include the matrix elements of the operators proportional to the total momentum of the 
QQ pair, the matrix element of Qi{^D2) and any other matrix element of octet operators 
with the exception of 08(^5"!), OsCSq) and OsC-Pi). In the power counting adopted in 
Q, which is described in [Q, all these matrix elements are suppressed by further powers 
of V and they can be neglected at this order of the expansion. For the same reason, the 
expression for r(35i —>■ l.h.) in does not include all the matrix elements of operators 
proportional to the total momentum of the QQ pair, the matrix elements of Qi{^Dj) 



-11- 



and Vi^^Si,"^ Di), and any other matrix element of octet operators with the exception of 
Os^Si), Osi^So) and OsePj). 

The P-wave decay widths at order v"^ are: 

repj - i.h.) = ^l^^^I^^{HePJ)\o^ePJ)lHePJ)) 

,2Im/8ePO^^(3p^)|0^(ip^)|^(3p^)^ + 2ImM^ 



^^ 2ImMi^.) ^^^3p^),g^(3^^)|^^3p^^^^ 2ImA( (ff (3p,)|0,(35O|i/(3p,)) 

k=l 

+ E ^.2 '^"^'-^^^''''^\ ff(^P2)|P.(^P2,^F,)|ff(3p,)), (2.30) 
j=l,8 

where J = 0, 1, 2. 

In Eq. ( |2.30| ), the first two matrix elements scale like v^, the remaining ones like v"^ . 



P-wave decay widths at order were computed in [10|, where the power counting of 



was used: they appear to contain only the first four terms of Eq. ( p. 30 ). It seems, however 



that also by adopting the power counting of Q at least the matrix elements of the operators 
'P8(^'S'i,^ -Di) and T'sacm should be added. 

3. Matching 

In this section, we calculate the order contributions to the imaginary parts of the 
matching coefficients that appear in Eqs. ( ^.28 )-( 2.30| ). The method consists in equating 



(matching) the imaginary parts of scattering amplitudes in QCD and NRQCD along the 
lines of 0. 

In the QCD part of the matching, the ingoing quark and the outgoing antiquark are 
represented by the Dirac spinors u{p) and v{p) respectively, whose explicit expressions are 



where Ep = y^T^ + M^, and ^ and r] are Pauli spinors. In the NRQCD part of the matching, 
the ingoing quark and the outgoing antiquark are represented by the Pauli spinors ^ and 
r] respectively. 
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We will match singlet, octet and singlet-octet transition operators at order a^; to this 
purpouse we will consider both the scattering amplitudes QQ QQ and QQ g — s- QQ, 
with no more than two gluons in the intermediate states. 

In the center of mass rest frame, the energy and momentum conservation imposes the 
following kinematical constraints on the scattering QQ — > QQ, 

\p\ = \k\, p + p' = 0, k + k' = 0, (3.2) 

and on the scattering QQ g — > QQ, 

Ep + Ep> + \q] = 2Ek, p + p' + q = 0, k + k' = 0, (3.3) 

where p, p' are the ingoing and k, k' the outgoing quark and antiquark momenta, while q 
is the momentum of the ingoing gluon, which is on mass shell. 

The matching does not rely on any specific power counting and can be performed order 
by order in 1/M [^]. We will perform the matching up to order which is the highest 

power in 1/M appearing in Eqs. (2.28)-( 2.3[l| ). In practice, we expand the QCD amplitude 



with respect to all external three-momenta. Note that, in the relativistic expansion, the 
gluon momentum \q\ is proportional to (three-momenta) ^/M. In the matching calculation, 
therefore, the gluon three-momentum appears with an extra 1/M suppression with respect 
to the quark and antiquark three-momenta. In the case of the QQ g QQ scattering, the 
expansion in the gluon momentum may develop infrared singularities, i.e. terms propor- 
tional to 1/1^. These terms cancel in the matching, as expected, having QCD and NRQCD 
the same infrared structure. For a detailed discussion see [^]. In the hadronic calculation, 
individual diagrams that contribute to the imaginary part of the QQ g — > QQ scattering 
amplitude containing interactions between the gluon in the initial state and gluon propa- 
gators develop also collinear singularities, i.e. terms proportional to 1/(1 ± cos9), 9 being 
the angle between the incoming gluon momentum and the momentum flowing in one of the 
gluon propagators put on shell to get the imaginary contribution. These singular terms 
cancel in the sum of all diagrams. Finally, we expect that, since the matching does not rely 
on a power counting and scattering amplitudes do not have a definite angular momentum, 
the matching will determine more coefficients than needed in Eqs. ( |2.28| )- (|2.30 ). 



3.1 QQ to light hadrons: singlet matching 

The matching of the QQ — > gg (qq) — > QQ amplitude is performed by equating the sum 
of the imaginary parts of the QCD diagrams shown in Fig. |l| (taken by cutting the gluon 
propagators or the light quark propagators according to 1/k'^ —2it i 5{k'^)6{k^)) to the 
sum of all the NRQCD diagrams of the type shown in Fig. |2|. The first two diagrams in 
Fig. H contain both a color singlet and a color octet part, coming from the decompositions 

2iV, 2iV, ^34^ 

t-t' ® tH' = ^1, 1, - —t^ , 

while the other five Feynman diagrams contribute only to the octet part. 
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Figure 1: QCD Fcynman diagrams describing the amplitude QQ QQ at order 



The calculation of the box diagrams in Fig. |^ gives the matching coefficients of the 
dimension 6, 8 and 10 singlet operators proportional to the relative momentum of the 
QQ pair, hsted in Eqs. (|A.16|) - (|A.19| ) and (|A.26| )- (|A.28| ). We quote the coefficients of the 
dimension 6 and dimension 8 operators in appendix They agree with those calculated 



in y]. We refer to |14| and references therein for an updated list of imaginary parts 
of matching coefficients of dimension 6 and 8 four-fermion operators; some of them are 
known at next-to-leading order. For dimension 10 operators we find 



Im hi{^D2) 
lmh'XSo) + lmh'(CSo) 
Imgi(3Po) 



2 ai 



15 2Nc ' 
68 vr Cp 
45 2Nc ' 
ttCf 



-7 ' 



2Nr 



(3.5) 
(3.6) 
(3.7) 
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Figure 2: Generic NRQCD four-fermion Feynman diagram. The empty box stands for one of 
the four-fermion vertices induced by the operators hsted in appendix 0, Eqs. ( A.16 )-( A.21 ) and 

(^^)-(|a3i]). 



Imffi(3p2,'F2) 



8 ttg vr Cf 

20 TT Cf 

21 2Nr 



(3.8) 
(3.9) 



The four-fermion operators to which the matching coefficients refer are Hsted in appendix 

The coefficients relevant for P-wave decay widths at order v''' are ( |3.7D and ( |3.8| ) . They 
were first computed in |jlO|, but our results disagree with the ones reported there. Note 
that while Eqs. (|]|) and (U) agree in the QED hmit with the results of the QED 
limit of the results in |jlO| is in disagreement both with [12| and [^]. 

Equation ( |3.6| ) agrees with the one found in Q . By matching the diagrams of Fig. |l] we 
cannot resolve luih'^i^So) and lmh'((^ Sq) separately. These coefficients multiply operators 
that contribute to the corrections of the S-wave decay widths. Eq. (^]^) contributes to 
the leading order decay width of the singlet state of the D multiplet, which for charmonium 
and bottomonium has not yet been observed; it agrees with the result of [Qq]. 



3.2 QQ to light hadrons: octet matching 

The calculation of the diagrams in Fig. || provides also the coefficients of dimension 6, 8 
and 10 color octet operators. Again, since we work in the center of mass rest frame, we 
cannot obtain the matching coefficients of the dimension 8 and 10 operators proportional 
to the center of mass momentum. 

The coefficients of the dimension 6 and 8 operators are quoted in appendix p| and agree 



with those obtained in Q and |17|. The coefficients of the dimension 10 operators are new 
results of this work. We find 



Imh'sCSi)+lmh'^{^Si) = 


29 2 1 2 

108 ' ^ 108 ' 
— a^vrn/ + —a'^irNc, 


(3.10) 
(3.11) 


Imhsi^Di) = 




(3.12) 


Im/i8(3D2) = 


— o^irNc, 
30 " ' 


(3.13) 
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ImhsCDs) = ^a^TTN,, (3.14) 

2 A^^ — 4 
lmhseD2) = Y^«s^^]v^> (3-15) 

fia /V2 _ A 

lmh'seSo) + lmhU'So) = ^«sVr^^, (3.16) 

Imgsi^Pi) = -^a^.TTN,, (3.17) 
_/V^ — 4 

ImgsePo) = -7«s vr^^, (3.18) 

Im58('P2) = -5"svr^^, (3.19) 

Im58(3p„3F2) = -|a^^|^. (3.20) 

The four-fermion operators to which the matching coefficients refer are hsted in appendix 
3.3 QQ g to light hadrons 

We show in figures ^§ the diagrams that contribute to the QQg — > QQ scattering amphtude 
with terms with color content (8) 1 or 1 (8) t". The imaginary part of QQg QQ is 
computed considering all possible cuts of the gluon propagators. Diagrams in Fig. ^ and 
Fig. develop collinear singularities, that however cancel when all possible cuts are taken 
into account. In these figures, the cuts are explicitly indicated. 

In the matching procedure, the QCD amplitude is equated to the sum of all NRQCD 
diagrams of the type shown in Fig. |^ These are all diagrams of NRQCD with an ingoing 
QQ pair and a gluon and an outgoing QQ pair. They can involve four-fermion operators 
and a gluon coupled to the quark or the antiquark line, but also four-fermion operators that 
couple to gluons. Four-fermion operators that induce octet to singlet transitions on the 
QQ pair may be one of the operators listed in Eqs. ( [A. 24 ) and ( [A. 25 ), but also one of the 



four-fermion operators involving only covariant derivatives, which, despite being usually 
denoted as singlet (or octet) operators, couple to the gluon field through the term —it°'gA°' 
in the covariant derivative and therefore have a singlet-octet component. 

The calculation of the imaginary part of the QQg — > QQ scattering amplitude allows 
to find the matching coefficients of dimension 8 and dimension 9 singlet-octet transition 
operators and dimension 8 operators proportional to the total momentum of the QQ pair. 



It also allows to fix the individual coefficients appearing in Eqs. ( |3.6D and ( 3.16| ). As 
discussed in section |2.3.2| , the basis of operators that we chose contains as independent 
operators Ti_^(^Pj,^ Si), with J = 0,1, 2. In appendix]^, we give for completeness also the 
matching coefficients computed with the other possible choice of independent operators, 
1/2 {Q'si'Si) - Q'ieSi)), 1/2 {Q'^eSu'Di) - Q'^{^Su^Di)) and TlfeSi,^P). This sec- 
ond set allows to establish the individual coefficients of the operators appearing in Eqs. 
( 3.10 ) and ( 3.11| ), but it is less useful for the discussion of the P-wave decay widths. 
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Figure 3: Box Diagrams: the gluon in the initial state interacts with a fermion leg. The other six 
diagrams, in which the two gluon propagators cross, have not been displayed. 



The matching coefHcients are: 



ImsLs^So,^ Si) = ~al7r + 



Imh[i^So) = ^ 



1 a^imf 
12 N, 



Imsi.8('Si,i5o) = 0, 



Cf 



1 



2 J- 2 

a„7r- " 



10 

9 "^''2iVe 48 

1 r, Cp 1 

5"^"2iVc ' 48 

10 2 N'^-A 
9 ' ANc, ■ 

2 2 iV2-4 



z, Cp 1 



(3.21) 
(3.22) 

(3.23) 

(3.24) 
(3.25) 
(3.26) 
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Figure 4: Box Diagrams: the gluon in the initial state interacts with a gluon propagator. The 
other four diagrams, in which the two gluon propagators cross, have not been displayed. 



Im4!)'(35i,3p) 

Im h'K^Si) 
lmh[{^Si,^D,) 
Im/i'/(35i,='Z)i) 



1 2 -4 
(T vr 

1 2 A^c'-4 

— cy TT 

12 ' 4Af2 ' 

1 2 

"12"^^ 4iV2 
1 2 A^c -4 

— (T TT 

1 2 ^^c'-4 

(T vr 

4^ 4iV2 ' 



(3.27) 
(3.28) 
(3.29) 
(3.30) 
(3.31) 



Imti.8(3Po,'5i) 



7 n 



3 r, Cp 

--at-K—^ + I h 

2 ' 2iVc V240 192 iV, 



Imt,3(^P„^5,)=(^ + ^^)«s^- 



1 



107 n 



1 



25 n/ 



(3.32) 
(3.33) 
(3.34) 



The four-fermion operators to which the matching coefficients refer are listed in appendix 
The total momentum of the ingoing QQ being different from 0, the matching calculation 
for QQg QQ also provides the coefficients for the operators defined in ( |A.22 ) and ( |A.23| ): 



Im /l cm 

I™ fi cm 



2 Cf 

4 " 2Nr 



1 

— I 

4 
0, 



(3.35) 
(3.36) 
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Imgiacm 
Im^lfecm 

Ini giccra 
Im^Sacm 

Imysccm 



0, 
0, 



1 2 C'f 

-— a^vr 

4 ' 2A^. 



24 



a^Trnf, 



0, 



1 2 ^^c'-4 

4 ' 4iVc 



(3.39) 
(3.40) 

(3.41) 

(3.42) 
(3.43) 

(3.44) 



We have checked the matching coefficients ( p.35[ )-( |3.44 ) by repeating the calculation of the 
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Figure 6: Vertex corrections; tlie gluon in tlie initial state couples to a fermion leg or to a three 
gluon vertex. 



diagrams in Fig. |T] up to order 1 /Af ^ in the general frame 

^ 1 



P = -^Q+Pr, 



P 



q—Pr, k' = -q — kr. 



k = 
k' 



1 ^ r 

-q + kr, 



2" 2 
Equations (|3l2ll)-( |3.44D are original results of this work. 



4. Poincare invariance constraints 

We can use Poincare symmetry to obtain independent checks on some of the matching 
coefficients derived in the previous sections. Here we outline the procedure, following the 



method of Ref . 

NRQCD is constructed by expanding (and matching) QCD in the non-relativistic limit. 
As a consequence, while translations and rotations are still explicit symmetries of the 
NRQCD action, the explicit invariance of the QCD action under boost is lost in the non- 
relativistic regime. However, the boost invariance of QCD manifests itself in a nonlinear 
realization, constraining the form of the NRQCD Hamiltonian. 

The constraints posed by Poincare invariance on the bilinear sector of the NRQCD La- 
grangian have been studied extensively in |13| and [|^. The computation of the matching 
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Figure 7: Vertex corrections: tire gluon in the initial state interacts with a gluon propagator. 



coefficients ( p.35 )-( |3.44| ) completes our knowledge of the imaginary part of the NRQCD 
Lagrangian at order including four-fermion operators proportional to the total mo- 

mentum of the QQ pair, which, due to their suppression in v have not been considered 
before. Therefore, we can now study the constraints induced by Poincare invariance in the 
four-fermion sector of the NRQCD Lagrangian. We adopt the method described in ||l^ by 
constructing the generators of time translation space translations P, rotations J and 
boosts K inside the effective theory and by imposing that the commutation relations of 
the Poincare algebra are respected. Since rotation and translation invariance are mani- 
festly maintained in NRQCD, the commutation relations involving only P and J are 
trivially satisfied while the commutation relations involving the boost generators K impose 
restrictions among the matching coefficients: 



[P\K^] 



(4.1) 
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Figure 8: Vacuum polarization: the gluon in the initial state interacts with a fermion leg. 






Figure 9: Generic NRQCD four-fermion Feynman diagrams involving an ingoing QQ pair and a 
gluon and an outgoing QQ pair. The black box with a gluon attached to it and the empty box 
stand respectively for one of the four-fermion-one-gluon vertices and for one of the four-fermion 
vertices induced by the operators listed in appendix Eqs. (A.16)-( A.31). The black dot with a 
gluon attached to it stands for one of the quark-gluon vertices induced by the bilinear part of the 
NRQCD Lagrangian given in Eq. (|A.l|). 



[J\K^] 
[K\K^] 



-iP\ 
—ie'^^^j''. 



(4.2) 
(4.3) 
(4.4) 



The construction of the generators proceeds in the following way: P and J can be 
obtained from the symmetric energy-momentum tensor [18, 19]: 



P 



iD]^ + xH-iD]x + 



1 

2 L" 



(4.5) 
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J = j d^x^^ {x^ [-iD^ +0 V' + X"^ X (^-i^) + X + X [n"x,B^, (4.6) 

where {1^,11^^), {x,ix^) and, in the = gauge, (Aj,!!^ = 5/:nrqcd/9(9o^")), are the 
pairs of canonical variables. The NRQCD Hamiltonian density /inrqcd can be obtained 
from a Legendre transformation of the Lagrangian density: 

-f^NRQCD = J C^^a^%RQCD 

1=1,8 

i=l,8 

The coefficient ci is equal to 1 at all orders in as, see |13] and |18|. 

A way to construct K is to write down the most general expression consistent with 
the NRQCD symmetries and to match it to the QCD boost generator, K = —tP + 

J d'^x— {x, /iqcd}- This procedure is analogous to the one followed in the construction 

of the NRQCD Lagrangian: new matching coefficients, typical of K, appear. The form of 
K in NRQCD is 

K = -tP+ I d^x]- {x, %rqcd} -f^l d'x^^&\ (4.8) 
■' 1=1 

This form is chosen in analogy to the QCD boost generator and satisfies ( [4.1| ) . ^^'^ contains 
all the possible operators with mass dimension M' that are vectors under rotation, are odd 
under parity and are invariant under C and T transformations. 

We now compute the imaginary, four-fermion part of the commutator ( [4. 21 ) at order 
To this aim, we need the bilinear NRQCD Hamiltonian at order 1/M, the four- 
fermion part of the NRQCD Hamiltonian at order the operator K^^^ and four- 
fermion operators in the boost generator, which first appear in ^W. The form of 

Kd) = i^t| X i-W)^ - ix^l X {-iD)x^ (4.9) 



and its coefficient ki were obtained in [|T8|, where it was shown that /ci = 1 to all orders in 
a^. In K'^^\ four-fermion operators like 

i?(4) = l^t| X {-iD)xxH 

appear. We do not give the detailed form of K^^'' since an explicit calculation shows that 
l/M^ j d^x[K^^\x),H] = 0{1/M^). 
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Using the canonical commutation relations we find for singlet operators at order 1/M^: 



X 



+ 



+ [d\^'^a\)x^a^i, - ^ta^-x5^(xVV)) (2Imgi,,^) 

-ie^'- (i,^o'"d^xxH - V'^XxW^V') Qlm/i(i5o) - Im/i. 



(4.10) 



and for octet operators: 



=^ld'x (d^{ijU-x)x^t-i^-i^^t-xd^ix^t-^)) Qlm/8(i5o) + 2Im58c 
+ [d^^H-a'x) X^iVV - V^tV^x9^(x^tVV)) Qlm/8(35i) + 2Im58,,^) 

+ (a^(V^^7^t'^x) X^^^V' - V^^T^x5*(X^tV») (2Im 586cm) 

-ie^''™ (^V^tV'V™xX^i> - V^t"x X'^i V'V^V) Qlm/8(^5o) - Im/8 

-ie^'^ (V^t^a^/^^t^a^™^ _ ^t^a^^^^t^a^Z^^ Qlm/8(=^5i) - Im 



cm 



0. 



(4.11) 



Equations (|4.10| ) and ( 4.11 ) imply that 



Im 5lccm 
Ini5'8ccm 
Imgiacm 
Im5(8acm 



Im/i(i5o), 
Im/8(^5o), 



4 
1 

'4 

-^Im/i(35i), 
1^ -.3, 



rIm/8("5: 



Imbibe 



Im58b 



(4.12) 



cm — ^) 
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Im/icm = ^Im/i(^S'o) 



---ImheSi), 
Im/8cm = ^Im/8(^5o), 



(4.13) 



Im/scm = 7lm/8(^5'i) 



Relations of the same form as Eqs. ( 4.12 ) and ( 4.15 ) hold also for the matching coefficients 



of the electromagnetic operators. Equations ( 4.12 ) and ( 4.13 ) imply that the knowledge 
of the imaginary part of matching coefficients of the dimension 6 operators completely 
determines the imaginary part of the coefficients of the operators defined in Eqs. ( A.22| ) 
and ( [A. 23 ), proportional to the total momentum of the QQ pair. The coefficients ( ^.351) - 
(3.44), obtained in the previous section, satisfy Eqs. ([4.12|) and (4.13). 



5. Summary and Outlook 

In the paper, we have calculated the hadronic inclusive quarkonium decay widths in 
NRQCD at order in the relativistic expansion and at order a^. The electromagnetic S- 



and P-wave decay widths have been previously calculated at order v'^ in |T^, 12 1. If we 
count as{M) ~ v"^, terms of order as(M)^?;^ are part of the next-to-next-to leading order 
(NNLO) corrections to the pseudoscalar S-wave decays and part of the NLO corrections to 
the vector S-wave and the P-wave hadronic decays. 



The results for the S-wave hadronic decay widths are given in Eqs. (|]2|) and ([23s| ) 



with the coefficients at order listed in appendix Let us first consider S-wave vec- 
tor decays. In the power counting of Q, those coefficients together with previous results, 
including contributions to the decay width coming from three-gluon decays and loop cor- 



rections [Q, |1^, |2^, 21, |2^, provide us with the full NLO expression of the hadronic 

inclusive decay widths, i.e. with the full expression up to order a'^v^, a^v^ and a^v"^ . In 

, , , ^ 21mdi^\^Si,^ P) 
the more conservative power counting adopted here, the octet terms AfS 

k=0,2 

x{H{^Si)\vll:l{^Si,^ P)\H{^Si)) need to be included. The matching coefficients d'i\^Si, 
^P) are however unknown. In the case of S-wave pseudoscalar decays, the largest uncertain- 
ties in the decay width come from the NNLO correction in to the matching coefficient 
Im/i(^S'o), from the NLO correction in to the coefficient lmgi{^So), and from the 
expression of Imd8(^S'o,^ Pi), which are all unknown. If we count as{M) ~ v"^, these are 
the only missing ingredient to complete the NNLO corrections to the pseudoscalar decay 
widths. Note that to complete the NNLO corrections to the pseudoscalar width, the NNLO 
expression of Im/i(^5o) and the NLO expression of Imgi(^S'o) would be necessary also in 
the power counting of 

The result for the P-wave hadronic decay width, calculated up to order as{M)'^v'^ , is 
given in Eq. ( p. 30 ) with the coefficients at order given in appendix In the case 
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of P-wave vector decays, the present calculation together with previous results, including 
contributions to the decay width coming from three-gluon decays and loop corrections 



1 17 , 24 ] , provides us with the full expression of the hadronic inclusive decay widths up to 
order afv^ and a'^v'^ . Explicitely we have: 
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i + ^n;) (/f(3p2)|Q^(35i,3Z)0|i?(^P2)) 

\ + Y2''f) (^('^2)|Q8('^l)|^('P2)) 



+- 



1 a^vr 



5 M6 



2 a^vr 



(^(^P2)|Q8(^I)2)|i^(^P2)) + ^^m'P2)m'D,)\H{'P2)) 



-^^{H{^P2)\VieP2:'F2)\H{^P2))-^^^{H('P2m{'P2,'F2)\H{^P2)). (5.3) 

A general source of concern is the proliferation of matrix elements with the increasing 
order of the expansion in v. Spin symmetry and vacuum saturation Q may help to reduce 
the number of independent matrix elements by relating different spin states and hadronic 
with electromagnetic matrix elements. The actual number of independent matrix elements 
depends on the power counting. 

In the power counting of only the first six matrix elements of Eqs. (^]^) and ( [5.3| ) 
and the first five of ( |5.2| ) contribute. In [0], it was assumed that only the first four matrix 
elements of Eqs. ( |5.ip and ( |5.3| ) and the first three of ( |5.2D contribute. 

The conservative power counting adopted here has been suggested in ||2^ to be appro- 
priate when Aqcd ^ mv'^. Under this condition, matrix elements are non-perturbative 
quantities and should be evaluated on the lattice. One can also take advantage of the fac- 
torization provided by potential NRQCD |26, 25, 11]. According to it, the matrix elements 
can be factorized into the product of the quarkonium wave function in the origin squared 
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(or derivatives of it) and few universal non-perturbative correlation function, eventually 
achieving a reduction in the number and a simplification of the non-perturbative operators 
needed. 

We also note that the convergence of the perturbative series of the matching coefficient 
is typically poor. For a discussion and references we refer for instance to |27]. 

Phenomenological applications of the expressions of the decay widths will therefore 
entail work in two complimentary directions: (1) improving the knowledge of the NRQCD 
matrix elements either by direct evaluation, for example by fitting the experimental data, by 
lattice calculations, and by models, or by exploiting the hierarchy of scales still entangled in 
NRQCD using EFTs of lower energy, like potential NRQCD; (2) improving the convergence 
of the perturbative series of the matching coefficients by resumming large contributions 
either related to large logarithms, or of the type discussed, for instance, in |28]. 
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A. Summary and definition of the NRQCD operators 

The two-fermion sector of the NRQCD Lagrangian relevant for the matching discussed in 
section ^ is: 

~ 1^ ° + 2M + 2M + 8M2 8M2 + 8M3 + 8M3 

32M4^ . L >y 32M4^ ^ ^' 16A/PJ^ 

+ c.c. , (A.l) 

where cr* are the Pauh matrices, iDq = ido — t"-gA(j, iD = iV + t"' gA"-, [Dx,E] = D x E — 
E X D, E^ = and B^ = {e^s = 1). We have not displayed terms of order 

or smaller and matching coefficients of 0(as) or smaller. The general structure of 
the four-fermion sector of the NRQCD Lagrangian is: 

n 

Here, we list the operators relevant for the matching performed in section ordered by 

< — > — > i — 
dimension. We use D = D — D . 
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For the octet operators defined in Eqs. (|A.2q) , (|A.2lD , ( |A.22D and (|A.29|) -( |A.31 



since 



the covariant derivative D does not commute with the color matrix t"" we need to specify 
the ordering between the two and verify that the resulting operator is gauge invariant. Let 
us consider, for example, the operator 08(^^1): 



and the three different orderings: 



X 



X 



Under the gauge transformation 



(1) 



(2) 



(3) 



Df^y X + '^^f'Dx, 
Dipy fx + iph^Dx- 



ij^ (1 + iw'^t'^) V', 
X^il + iu^'^nX, 

9 



UJ 



( A. 4 ), ( A. 5 ) and ( A. 6 ) transform respectively as: 



i'^Dt'^X 



X 



X 



(2) 



(3) 



(2) 



fabc^c ^\Dt^x - r''^duj^Tph^x 



jabc^c ^tDfb^ 



(3) 



Only the last ordering leads to a gauge invariant definition of OsCPi)'- 



6 0sePi)=r'"'u;' 
Therefore, we define 



i^^Dt'x 



(3) 



(3) 



+ 



(3) 



X^Dt^'^P 



(3) 



^t^^a^ = - (i5v^) ^ ex + i^U^'Dx. 



(A.3) 

(A.4) 
(A.5) 
(A.6) 



(A.7) 



(A.8) 
(A.9) 
(A.IO) 



= 0. 
(A.ll) 

(A.12) 



Generalizing to operators containing more than one covariant derivative, we define 



t 



V ... D '"f'x = (-1)" [D'' ■ ■ ■ D'^^j f'x 



(A.13) 



X- 



The singlet-octet transition operators are denoted by Oi.si^^^^Lj,'^^'^^ L'j,) or Os-i 
(25+iLj,2^'+iL'^,)- In the first case the first set of quantum numbers refers to the QQ pair 
in a color singlet state, the second to the QQ pair in the octet state, while in the second case 
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the first set of quantum numbers refers to tlie QQ pair in a color octet state and the second 
one to the QQ pair in the singlet state. In some cases, it has been found convenient to 
introduce singlet-octet transition operators that annihilate (create) states containing a QQ 
pair and a gluon in which the total angular momentum J of the the quark-antiquark pair 
has not a definite value (a definite value could be attributed by further decomposing these 
operators in irreducible spherical tensors). This is the case of the operators Ti^i^Si,^ P) 
and T^^\{^Si,^ P). In these cases, we cannot use the quantum number J and we have to 
denote the state just by the orbital angular momentum and spin quantum numbers. 

The symbols A^'^B^'^ and S^'^'^^'^ A^\ used in the definitions of some four-fermion opera- 
tors denote symmetric and traceless two and three indices tensors, according to 

A^'B^^ = ^ - —A ■ B, (A.14) 
2 3 

^ (A.15) 
For some details on the decomposition of Cartesian tensors in terms of irreducible spherical 
tensors see |12]. 



Operators of dimension 6 



OsC^o) =^h'^xx^t''^ij, 



(A.16) 



(A.17) 



Operators of dimension 8 



Vl{^Sl)=\^l,^ax■X^o(^-'-"D^ V + H.c. , (A.18) 



OieP2) (-^^V)) xx^ (-^'d(v^)) ^. 



(A.19) 
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OsCPi) =\^^ [~'{d X a) ex ■ X^ [-'{d x 
08('P2) =V'^ (^-^'d(V^)) exX^ (^_l^(vi)^ 

01 cm = {~\^) ^ ■ ^ (^^^^) + ' 

Oscm = -V'^ (-1"^) ■ ^ (x^i'^V-) + H.C. . 

Vlacm = V* (V^V^X) (XV^V) , 
^l&cm = V • {^J^ox) V • (x^'?V') , 
Vlccra = V (V't^) • V (xV) , 

T'sacm = Dl, [^J^ea^x) D'L (x^tV^) , 

Vsbcra = Dab ' (i^^t^^x) Dac ' {x^t'a^) , 
ccm — 

5i.8('5o,3 5i) =^V^5^ • axxH + H.c. , 
^i.sl^'Si,! So) =l^^9Bx ■ X^a^ + H.c. . 



Operators of dimension 9 

T,.sCSo,' Pi) =^4^^xxH^ -gE + gE- "D)nj + H.c. , 
JPg-sC^o,' Pi) =\r'"''4^h''xX^t\'D ■ gE- + gE'^ ■ "d)^ + H.c. , 
P8-8('5o/ Pi) =ld^'"'i^h\x^t'CD ■ gE- + gE' ■ + H.c. , 
V^lCS,,^ P) =^d-'^^^l;H-ax ■ X^a (d ■ gE^ + gE^ ■ ^)t'^ + H.c. , 
V^li^Si,^ P) =]^(F^''i,h''a\x^(r^ (D^'gE^^^ + gE^^^ ^^)t''i^ + H.c. , 

^&{^So) =lr'"'i^\D . EffxX^e^ + H.c. , 
ri.8(^Pi/ -So) =\i^''9Ex ■ X^^Di^ + H.c. , 

r/°) P) J-^^&x ■ X^B{D .gE + gE- ^)ip + H.c. , (A.25) 
tI}^{^Si,^ P) =^iP^ax ■ X^a x {-"d x gE - gE x "d)^ + H.c. , 
Ti5'(35i,3 P) =^^j^ax ■ X^a x (d x gE - gE x "d)^, + H.c. , 
T^h^Si? P) =\^h''ux ■ X^B x(d X gE<^ - gE'^ x "d)^ + H.c. , 
Tl^^eSi,^ P) =]^'il;^a'xX^a^CD^'g&) + gE^'lD^^)^; + H.c. , 

ri.8(3po,' Si) =^v^ p • <?) xx^<? • gE^ + H.c. , 

ri.8('Pi,=' Si) =]ijj^ ("d xa)x- X^d X gE^P + H.c. , 



Ti.s{'P2,^Si) =^V^ ('d(VJ)) xxV^*5^'V + H.c. 
Operators of dimension 10 

Q'li'So) (-'-d)\x^ (-'^d)\, 



Q"i(%)=^V'^(-|'^) xxV + H.c, 



Q'ii'Si)=i^^\^--Dj ax-X^[-^Dj aV, 
Qii'Si) =^V^ (^-^"d y ax • X^'?V' + H.c. , 

Q'/(=^^i,=^i^i) =^V'^ {-l^y (-0'^(^^^V^xxV^ + H.c.. 



(A.26) 
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=^^^ (-^^^V^'^) (-^^)^^' (~^^(V)) V + H.c, 



Ql(3£)2) =^^t |^giim^(i^Z)^m ^ 1 ^ij7^(m^0(jm^ ^ (A.28) 



(A.27) 



Q^e^i) =V'^ (-l^)^i'^X-X^ (-^^^atV, 



Q8r5i)=2V'M-2^ ) <?i"x-X^<?i"V' + H.c., 



2 / . \ 2 



'o^'^D^^a' i"x X'^f^-'i" + H.c. . 



(A.29) 
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(A.30) 



B. Summary of matching coefficients 

In the following, we list all the imaginary parts of the matching coefficients of the four- 
fermion operators up to dimension 10, calculated at 0{a^) in the strong coupling constant 
in section |^. 

In the presentation of the results we give for completeness also the matching coef- 
ficients obtained by using a basis of operators that includes 1/2 [Q.'^{^ Si) -Q!^{^Si)), 
1/2(Q^(35i,3Di) -Q!l{^Si,^Di)) and T^}l' Si,"" P) instead of TlsI^Pq," Si), ^i-sl'^i, 
'7i_8(^-P2 5^ 5'i). It is understood that when this basis is used, the coefficients \mti_^(^ Pj , 
with J = 0, 1, 2 are set to 0. Viceversa if our basis contains the operators Ti.^i^Pj,^ Si), 
with J = 0, 1,2, the coefficients Irnh'^i^Si) -lmh'^{^Si), Irnh'^i^Si,^ Di) -Imh'^{^Si,^ Di) 
and Im4!l'(^5'i,^P) are set to 0. 



Operator of dim. 6 Matching coefficient Im (Value) 

OiCSo) Imfii'So) a^^^ 
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Oi('Si) 



Im/i(35i) 







TV^ — 4 



Operator of dim. 8 Matching coefficient Im (Value) 



4 2 Ci. 
3 ' 2Nc 



OiCPi) Im/i(iPi) 



2N, 

Oi('Pi) Im/i(3pi) 



'5i-8(^5'o,^ 5*1) 
'5i-8('^5'i/ So) 



Imsi_8(^S'o,^S'i) 
Imsi_8(^5'i/ 5o) 



a^TT ( 1 



4iVc V3 







Irngs^Si) 
Irngsi'Su^Di) 



4 2 - 
3 ' ANr 



9 ' ^ 



6 ' ^ 



i, 
i, 
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Im/8(^Po) 
Im/8(=^P2) 



12 

SagTr — — — 


-Qf^TT — £ 

5 ' ANr 



i, 



1 2 C'f 



Clem Im/icm 4"''^2A^ 

(-^8 cm im/scm T"s^' 



1 



4 4A^c 



"Pla cm Ini ffla cm 

'Plfecm Imglfecm 

T> T 1 2 C'f 



4 " 2A^c 

1 ,2„ 



^8acm Imgsacm -—a^TTUf 

^86 cm Im58fecm 

1 o iV2_4 



'P8ccm Im^Sccm -T^s^ 



4iVc 



Operator of dim. 9 Matching coefficient Im (Value) 

T^^'i'Sr,' P) ImiSr Si,^ P) -r"^^^ 

r^'i'S,:^ P) Imt^rSi,^ P) ^^al.^^ + logV - iag^vrg 



T,3(^Po,^ ^0 Imt,srPo,^ 50 --ag^vr^ + I 240 + 192 it ' 



Ti.8(^Pi/5i) lmti.s{-'Pi;'Si) — + 



V72 576iVcy 
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10 576 iVc. 



Operator of dim. 10 Matching coefficient Im (Value) 

1 /v2 _ 4 



Vii'Pi) 

Vl(^P2) 

rieP2,^F2) 



Imr7i(iPi) 

Im^iCPo) 
Im5i(3pi) 

Im5i(3p2) 

Im 51 (3^2,3^2) 





2Nr 





89 C'f 

— a„TT 

5 ' 2Nc 

20 2 

— a„7r 

21 ' 2iV^ 



Im/iiCL»2) 

Im/ii(3Di) 
lmhi{^D2) 



15 



2 

^"="2iV^ 






Imhg(}So) 



2 /v2 _ 4 



5 ^ 4iV, 



c 



+ i^;,^(35^,3^^) + i^^.(3^^^3^^) |«^vrn, + la^TTAT, 



7V^ — 4 



47Vc 

VsCPi) IragsCPi) 

VsCP2,^F2) Im58(3p2,'F2) __a2^^__ 



9 /v^ _ 4. 

Q8('^2) lnihsCD2) —aln- ' 



15 ' AN, 



c 



21 



References 



[1] N. Brambilla et al, "Heavy quarkonium physics," CERN-2005-005, (CERN, Geneva, 2005) 
[arXiv:hep-ph/0412158]. See also the web page of the International Quarkonium Working 
Group : http : / / www. qwg . t o . infn . it . 

[2] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Rev. Mod. Phys. 77, 1423 (2005) 
[arXiv:hep-ph/0410047]. 

[3] W. E. Caswell and G. P. Lepage, Phys. Lett. B 167, 437 (1986). 



-38- 



[4] G. T. Bodwin, E. Braatcn and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. 

D 55, 5853 (1997)] [arXiv:hep-ph/9407339]. 

[5] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 
[6] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 
[7] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592, 1 (2004). 
[8] A. Vairo, AIP Conf. Proc. 756, 101 (2005) [arXiv:hep-ph/0412331]. 
[9] G. T. Bodwin and A. Petrelh , Phys. Rev. D 66, 094011 (2002) [arXiv:hep-ph/0205210]. 
[10] H. W. Huang, H. M. Hu and X. F. Zhang, Phys. Rev. D 56, 5816 (1997). 

[11] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Lett. B 580, 60 (2004) 

[arXiv:hcp-ph/0307159]. 

[12] N. Brambilla, E. Mereghotti and A. Vairo, JHEP 0608, 039 (2006) [arXiv:hep-ph/0604190]. 

[13] A. V. Manohar, Phys. Rev. D 56, 230 (1997) [arXiv:hep-ph/9701294]. 

[14] A. Vairo, Mod. Phys. Lett. A 19, 253 (2004) [arXiv:hep-ph/0311303]. 

[15] J. P. Ma and Q. Wang, Phys. Lett. B 537, 233 (2002) [arXiv:hep-ph/0203082]. 

[16] V. A. Novikov, L. B. Okun, M. A. Shifman, A. L Vainshtein, M. B. Voloshin and 
V. L Zakharov, Phys. Rcpt. 41, 1 (1978). 

[17] A. Petrelh, M. Cacciari, M. Greco, F. Maltoni and M. L. Mangano, Nucl. Phys. B 514, 245 
(1998) [arXiv:hep-ph/9707223]. 

[18] N. Brambilla, D. Gromes and A. Vairo, Phys. Lett. B 576, 314 (2003) 

[arXiv:hep-ph/0306107]; N. Brambilla, D. Gromes and A. Vairo, Phys. Rev. D 64, 076010 
(2001) [arXiv:hep-ph/0104068]. 

[19] A. Vairo, Nucl. Phys. Proc. Suppl. 133, 196 (2004) [arXiv:hep-ph/0310064]. 

[20] P. B. Mackenzie and G. P. Lepage, Phys. Rev. Lett. 47, 1244 (1981). 

[21] H. W. Huang and K. T. Chao, Phys. Rev. D 54, 3065 (1996) [Erratum-ibid. D 56, 7472 
(1997 ERRAT,D60,079901.1999)] [arXiv:hep-ph/9601283]. 

[22] M. Gremm and A. Kapustin, Phys. Lett. B 407, 323 (1997) [arXiv:hep-ph/9701353]. 

[23] J. Campbell, F. Maltoni and F. Tramontano, Phys. Rev. Lett. 98, 252002 (2007) 

[arXiv:hep-ph/0703113]. 

[24] K. Hagiwara, C. B. Kim and T. Yoshino, Nucl. Phys. B 177, 461 (1981); R. Barbieri, 
E. d^Emilio, G. Curci and E. Remiddi, Nucl. Phys. B 154, 535 (1979). 

[25] N. Brambilla, D. Eiras, A. Pineda, J. Soto and A. Vairo, Phys. Rev. Lett. 88, 012003 (2002) 
[arXiv:hep-ph/0109130]; N. Brambilla, D. Eiras, A. Pineda, J. Soto and A. Vairo, Phys. 
Rev. D 67, 034018 (2003) [arXiv:hep-ph/0208019]. 

[26] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566, 275 (2000) 
[arXiv:hep-ph/9907240]. 

[27] A. Vairo, arXiv:hep-ph/0212271. 

[28] G. T. Bodwin and Y. Q. Chen, Phys. Rev. D 64 (2001) 114008 [arXiv:hep-ph/0106095]. 



-39- 



